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We construct a tensor network that delivers an unnormalized quantum state whose 
coefficients are the solutions to a given instance of 3SAT, an NP-complete problem. The 
tensor network contraction that corresponds to the norm of the state counts the number 
of solutions to the instance. It follows that exact contractions of this tensor network are 
in the #P-complete computational complexity class, thus believed to be a hard task. 
Furthermore, we show that for a 3SAT instance with n bits, it is enough to perform 
a polynomial number of contractions of the tensor network structure associated to the 
computation of local observables to obtain one of the explicit solutions to the problem, if 
any. Physical realization of a state described by a generic tensor network is equivalent to 
finding the satisfying assignment of a 3SAT instance and, consequently, this experimental 
task is expected to be hard. 



1 Introduction 

The rational to develop a quantum computer follows from the fact that simulating quantum 
systems in a classical computer is a hard task pQ. Nonetheless, an enormous progress has 
been achieved to improve on the classical simulations of relevant quantum problems using 
Monte Carlo and Tensor Networks techniques. In particular, networks of tensors can be used 
to represent the coefficients of the wave function in a manner that faithfully represents the 
entanglement of the quantum state. The most relevant tensor networks geometries correspond 
to one dimensional Matrix Product States (MPS [2j [3]), to higher dimensional Projected 
Entangled Pair States (PEPS [4j [5] ) and to Multiscale Renormalization Group Ansatz (MERA 
[6]). The limits for the efficiency of the tensor networks technology depends on the amount 
of entanglement pervading the system. 

The main difficulty encountered in the application of tensor networks to simulate quantum 
states corresponds to computing observables, that is, to perform the contraction of all the 
indices in the tensor structure. Let us illustrate this point using the computation of the 
properties of the ground state of a given Hamiltonian. The tensor network approach needs 
two steps. First, it is necessary to find the individual tensors in the tensor network that 
give the approximation to the ground state. This can be done using different strategics. A 
possibility is to use local Euclidean evolution, based on the idea of using the Trotter expansion 
for the total Hamiltonian divided in local terms (see for instance ITEBD for infinite spin 
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chains [7]). Second, once the tensor network has been constructed, a contraction of the tensor 
network is needed to perform the computation of an observable. This contraction can be 
shown to be efficient in one dimension but becomes unpractical in higher dimensions due to 
the typical area law scaling of the entropy found in any quantum system governed by local 
interactions. Approximate contractions of tensor networks do produce acceptable results, 
though the problem is worse and worse as the connectivity of the system increases. 

In order to assess rigorously the complexity class associated to the problem of contracting 
arbitrary tensor networks, we here present an exact representation of the solutions of a 3SAT 
problem as an explicit tensor network. As we shall see, the tensor network that delivers 
the solutions to a 3SAT instance can be constructed analytically without the need of any 
minimization or training process. It is then possible to relate the problem of contraction 
of tensor networks to complexity classes associated to the 3SAT problem. As a matter of 
fact, the 3SAT problem has been addressed previously using quantum strategies [8]. In this 
context, the problem was analyzed using MPS in Ref. [9j 110) as an approximation technique 
to analyze the typical gap in an adiabatic computational approach to NP-complete problems. 
Big instances of the problem can be solved, however, using classical statistichal methods 

mug. 

3SAT is an NP-complete problem p~3j [14] , the complexity class of decision problems ver- 
ifyable in polynomial time. Finding an explicit solution to an NP problem defines the FNP 
complexity class, while counting the number of solutions to a polynomially computable func- 
tion defines the #P complexity class. For instance, counting the solutions of an instance of 
3SAT is in #P, thus believed to be a hard task. We show that the contraction of tensor 
networks for the 3SAT problem is a #P problem, and obtaining an explicit solution to the 
problem and reducing the bond dimension of the tensor network, are both in FNP. 

This paper is organized as follows: in Section [2] we present an explicit construction of 
the tensor network, used in Section [3] to analyze the complexity of the operations related to 
tensor network manipulation. In Section |4] we show that reducing the virtual index of a tensor 
network is a hard problem. Finally, we draw the conclusions in Section [5] 

2 Explicit construction of a tensor network for 3SAT 

Let us start by recalling the definition of a 3SAT problem. An instance of 3S AT corresponds to 
a decision problem, namely, determining whether there is an assignment of a set of n bits that 
can take values Xi = 0, 1, where i = 1, . . . , n such that a number of constraints are satisfied. 
These constraints are defined by the set of m clauses (C a ,a = 1, . . . , m), where each clause 
involves three bits and rejects one of the eight possible assignments. To be more explicit, 
let us consider a clause C a involving bits i, j and k and take the case where the assignment 
(yi,yj,Vk)a is excluded. This amounts to impose the boolean restriction yi A yj A y*,. The 
3SAT problem is defined as the problem of deciding whether there exists a string of bits 
satisfying all the clauses. As said, the 3SAT problem is NP-complete. 

Given a 3SAT instance, we construct a graph adapted to this specific instance under 
consideration. The graph contains two kind of structures: one is associated to bits and another 
one to clauses. Each bit carries an open index related to Xi and indices linking the bit to the 
clauses related to it. Clauses have no open index, they simply connect to the bits they relate. 
An example of such a graph is presented in Fig. [T] This kind of graph naturally represents 
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the connectivity of the 3SAT instance and is common to many approaches to 3SAT as e.g. 
the belief propagation strategy pTJ [12] . It is important to notice that not all instances of the 
3SAT problem are hard. Considering a problem over n variables and m clauses, and defining 
a = m/n, the problem is believed to be hard only in the critical range 3.42 < a c < 4.506 
[TT] . This constraint restricts the connectivity of tensor networks descriving hard instances of 
3SAT, which are represented by tensor networks only slightly connected. In our construction 
we only require m + n tensors and 3m links. Notice that connectivity, throught the planarity 
of the network, identifies as well hard instances of the Ising model in higher dimensions |15j . 




(circles) and clause tensors (squares). Each clause tensor is connected to 3 bit tensors 

through the ancillary bonds a, and 7. Bit tensors can be connected to an arbitrary number of 
clause tensors, and possess a physical index related to the variable x^. The contraction of this 
tensor network provides a solution to the 3SAT instance. 



Our next step is to construct an unnormalized quantum state based on this structure. 
The idea is to associate a tensor to each element of the connectivity graph presented above. 
Each bit is now associated to a tensor Qa^ fc .., Q ,. fc that describes qubit k. The index can 
take values or 1, that is the two possible states of the qubit. The tensor for this qubit 
also carries ancillary indices that we have represented with Greek letter a's connecting to rj, 
different clauses. Each qubit will connect to a different number of clauses, so the total number 
of links rfc will depend on the qubit k we are considering. Every ancillary index runs up to a 
maximum value \, on = 1, . . . , X- Then, we further introduce a tensor for each clause C a , that 
we represent as , where a = 1, . . . , m. Note that a, (3 and 7 are ancillary indices that 
connect to three different qubits. The contraction of all ancillary indices define our global 
unnormalized state 

\i>{Q,C))= t h >-^\ii,...,i n ) (1) 

i\ ,...,i n 
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where the coefficients t* u '"' t ' 



correspond to the contraction 



t 



(U, ■ ■ ■ , in\i>(Q, CO) = Tr . . . QNi-cW . . . <?M 



(2) 



where the symbol Tr represents the full contraction of all the ancillary indices that have been 
omitted in the above expression for the sake of simplicity. 

The problem we are addressing is how to assign the tensors Q^ lk and such that the 
coefficients of the tensor t are equal to 1 for the solutions to the 3SAT instance represented 
by the graph and otherwise. That is, 



Then, it follows that the unnormalized state \ip) corresponds to the equally weighted super- 
position of all the solutions to the 3SAT instance. 

It turns out that there is a very simple tensor network, with only % = 2, that carries all 
the solutions to a given 3SAT instance. The basic idea comes from the fact that such a tensor 
can be thought to be the ground state of a Hamiltonian made by the addition of 3-body 
Hamiltonians, each one associated to one of the clauses that define the instance, 



a=l 

The 3-body Hamiltonian associated to the clause C a involving qubits k is constructed as 
a projector on the wrong assignment to be rejected (jji,yj,yk)a, 



That is, for a given clause C a , its corresponding Hamiltonian H a penalizes with one unit of 
energy the wrong assignment to the clause and delivers a zero value for the rest of the allowed 
assignments. The relevant point is that this Hamiltonian is diagonal in the computational 
basis, and all the pieces do commute [H a , Hf,] =0, Va,6 = l,...,m, that is, the Hamiltonian 
is frustration free, and local minimization provides a global minimum. Under some circum- 
stances tensor networks i.e. PEPS have been identified as ground states of frustration- free 
Hamiltonians [16 . Therefore, we can look for a global tensor construction by just guar- 
anteeing the correct local tensors clause by clause. Hence, we obtain our first and central 
result: 

Result 1: The solutions to an instance of the 3 SAT problem can be represented explicitly as 
an unnormalized quantum state which is encoded in a tensor network with x — 2. 

The precise construction of the tensors that encode the solution of the 3SAT instance 
corresponds to 





rn 




(4) 



H a = \yt\vf,y^){yt\yf,yt ] \. 



(5) 





and 



C 




{ 



1, K A 7) = {VhVj,Vk) 
0, otherwise 



a 



(7) 
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All qubit physical indices ik are locked to the ancillary ones ati, and the values of the qubit are 
neutral, taking the same value 1 whether i = 1 or 0. The tensor for the clauses reads the qubit 
value through the ancillary indices and simply provide a 1 for all the allowed assignments but 
a for the rejected one. This structure is easily seen to minimize every H a and, thus, provides 
the global ground state which is the even superposition of all the global valid solutions to the 
set of clauses. That is, the contraction of tensor network we have constructed does fulfill Eq. 

El 

From our construction it is clear that there is a degeneracy of the tensor networks delivering 
the same state It is possible to multiply any coefficient of e.g. a qubit tensor by an 
arbitrary number, while dividing it in the associated clauses. It is also possible to assign 
different lockings between the variable tensors and the clause tensors so to produce the same 
solution. This freedom is reminiscent of the so-called gauge symmetry also present in other 
tensor network geometries [5]. Local changes in the tensor values produce a continuum of 
tensor network constructions leading to the same solution. Our choice in Eqs. El [7] can be 
understood as a canonical form for this tensor structure. 

The assignment in Eqs. El [7] is also found when a numerical Euclidean evolution is per- 
formed to the flat distribution (|0) + 1))®". The Euclidean minimization of the Hamiltonian 
H produces the equally flat distribution of solutions. In some sense, this is surprising. It 
must be made clear that if a solution to the 3SAT instance exists, a simpler network with 
the qubit indices being locked to their solution values will deliver the correct solution. From 
an initial random assignment to the qubits, we recover as well all the solutions with different 
weights. In the case of a single satisfying assignment, a tensor of x = 1 (a product state) 
also provides the solution. Yet, the Euclidean evolution keeps restoring symmetry on qubit 
assignments. As a consequence, the x = 2 tensor network we have presented is obtained as 
an infrared fixed point of an Euclidean evolution minimization algorithm. 

To further analyze the performance of the Euclidean evolution, we perform a numerical 
simulation of the tensor network subject to the evolution operator. In Fig|2] we show the 
overlap between the final state encoding the solution and the evolved state along the imaginary 
time evolution. We fix St — 0.001, resulting in n — 1000 steps to provide the exact answer. 
The commutation between different H a ensures that every step of the computation actually 
projects into the correct answer, identifying monotonically the right assignments already in 
the first steps of the evolution. This is reflected on an identical evolution for instances with 
the same number of variables and solutions. 

3 Complexity class for the contraction of tensor networks 

Note that the tensor network we have constructed is not normalized. Indeed, if the 3SAT 
instance associated to this tensor has no solutions, all the contractions leading to the compu- 
tation of each t n '■■■>*« will produce a zero value. Such a state would have zero norm. In the 
case that the 3-SAT instance has p satisfying assignments we state the following: 

Result 2: Given the above tensor network that encodes the solutions to a 3SAT instance 
in a state \ip(Q,C)), the norm {ip(Q,C)\ip(Q,C)) =p reads the number of solutions p to the 
3SAT instance. Hence, computing the norm of the tensor network explicitly corresponds to a 
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1000 



Euclidean steps 



Fig. 2. For instances of 3SAT with n = 12, 14, 16 variables we plot the overlap of the evolved state 
\ipt) at time t with the final state \iI>f) that encodes the correct solution as described by Eqs.[6][7] 
along the Euclidean evolution. We set the imaginary time step St = 0.001. The initial overlap 
and the evolution depends only on the number of variables and the number of solutions to the 
instance. 



#P- complete problem. 

This result is derived from the tensor network construction we have presented. It shows 
that a contraction of the tensor network counts solutions of a 3SAT instance, which defines 
the #3SAT problem. Consequently, the problem of contraction of arbitrary tensor networks 
falls into the #P-complete complexity class. This result was derived in Ref. [17] in the contex 
of contracting PEPS with bond dimension \ = 2- 

Computing the number of solutions can be also reduced to the simpler case of contracting 
\ip(Q,C)) against a x — 1 tensor network representing the unnormalized state \iptest) = 



(|0) + \\))® n . This contraction computes (V>t«ri#(Q, C)) = Tr Q [1] . . . Q [n] C^ . . . 
with Q [k] = Q[ fc l° + QM 1 . That is, the simplest non-trivial contraction of a tensor network 
with x — 2 with another with x = 1 suffices to compute p. Note that this problem can also 
be seen as finding the overlap of the 3SAT state \ip(Q, C)) with the neutral product state 
latest) = J2 X \ x )- 

Let us emphasize further that the tensor network we have designed only uses tensors with 
X = 2 rank. This shows indirectly the enormous representation power of tensor networks that 
use small rank individual tensors. Important resources for quantum information processing 
can be constructed as well with tensor networks with x = 2, such as the cluster state [IB] or the 
toric code state [19] . It was known that x = 2 tensor network provide good approximations to 
the ground states of relevant quantum interactions. PEPS with small rank for x were shown 
to support area law scaling of entanglement and to be able to describe a scale invariant theory 
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where all correlators decay polynomially [19 . 

In our case, the situation is more radical due to the fact that the connectivity is non-local. 
Then, the maximum entanglement support offered by the tensor network can obey a volume 
law. More precisely, we can take a 3SAT instance and produce a bipartition of qubits into 
two sets A and B. Consider the case where all qubits in A are involved in a clause that 
further relates to qubits in B. Our tensor network structure will then have at least n/2 
ancillary indices connecting A and B. Thus, the amount of entropy which can be supported 
by the tensor network is Sa ~ log 2 ~ n /2, where Sa is the von Neumann entropy of the 
reduced density matrix for the subsystem A. This is the largest possible entropy. Thus, our 
tensor networks can support volume scaling for the entropy. Numerical evidence support this 
necessity [20} 1ST]. 

Let us now turn to the problem of finding an explicit solution to a 3SAT instance, which 
is a FNP problem. This can be achieved using contractions of the tensor network we have 
constructed. From the fact that the 3SAT problem is self- reducible [14], n contractions of 
tensor networks deliver an actual solution, which leads to the following: 

Corollary: An explicit solution to a 3SAT instance, encoded in a state \tp(Q,C)}, can be 
obtained by computing a polynomial number of successive expectation values on that state. 

We now turn to present a method to extract one of the possible solutions, if any, of a 3SAT 
instance tensor network. The basic idea is to perform a series of n contractions as follows. 
We first take any qubit 1 and compute 

_ MQ,C)|zi|V>(Q,C)) 
[1> (iP(Q,C)\iP(Q,C)) 

where z\ = + I), an operator that reads 1 or for the elements of the basis of a 

qubit Hilbert space. This operation is only defined if the norm of the state is different from 
zero, that is, if there is at least one satisfying assignment solving the 3SAT instance. The 
computation of this expectation value amount to a new contraction of the tensor network 
that needs exactly the same amount of computational cost as finding the norm of the state. 
The result of the contraction can be any number between or 1. If the result is either one of 
the extremes, namely or 1, that means that the qubit must take such a value. The classical 
bit associated to this qubit is completely fixed. In contrast, the case where the solution is not 
one of the two extreme values corresponds to having more than one solution to the instance. 
In such a case, the state \ip{Q, C)) is a superposition of various solutions. The qubit which is 
being analyzed cannot commit itself since it takes different values for those different solutions. 
Consequently, we can take any of the two possibilities and construct a new tensor network 
where Qal', lk .,a T is fixed either to or 1, as we prefer. Both values will lead to a valid solution. 
We can now proceed to a second qubit and repeat the operation. In this way, n contractions 
of tensor networks will suffice to deliver one of the possible solutions to the 3SAT instance. 

It should be noted that our approach requires a polynomial number of calls to #P contrac- 
tions. As a matter of fact, however, each call should only need to decide whether a solution 
exists, i.e. it would suffice to know if a certain assignment can be satisfied at each step. Thus, 
each call would only amount to solving a NP problem. 



(8) 



8 An exact tensor network for the 3SAT problem 



4 Reducing the rank of tensor networks 

A problem common to all the available techniques to make a contraction of a tensor 
network corresponds to finding the best approximation of a tensor network with large % using 
another one with a small x- For instance, in two dimensions, the contraction of a tensor 
network accumulates open ancillary indices as dictated by the area law of entanglement. The 
contraction of a tensor network needs a larger and larger effective rank \- To solve this 
problem, different reductions can be tried. It is, then, necessary to approximate large rank 
tensor networks by small rank ones. 

The problem of reducing tensor networks is simpler in one dimension. For MPS, the 
optimization problem of finding an appropriate reduction of the ancillary indices requires 
polynomial resources [22], though the global optimization problem can be extremely hard 
even in ID [23]. For the general of non-local tensor networks this task is out of reach. This 
can be expressed as an evident statement. 

Result 3: The reduction of a tensor network \ip(Q,C)) that encodes a 3SAT instance to a 
\ = 1 tensor network is equivalent to finding the solution to that instance. 

This limitation can be deduced from our above results. Let us illustrate this point in 
the following way. As mentioned above, given an instance of 3SAT with only one satisfying 
assignment, we know of at least two tensor networks that produce the correct identical state. 
One of them is the tensor construction \tp(Q,C)) with \ — 2 we have presented and a second 
one is a tensor network with x = 1 corresponding to the product state. Both states are 
identical. Furthermore, we know a priori the first construction since it is neutral in all the 
tensor assignments. Then, the problem of finding the best product state approximation to 
the x = 2 tensor network, that is finding 

max J(V> x =i|Vx=2)| (9) 

is equivalent to finding the solution to the instance. But, it was shown that finding the 
solution to the 3SAT instance is in FNP. Thus, finding the optimal more economic tensor 
network construction of a state given a more costly one is also in FNP. 

Let us briefly discuss the relation of this result with practical applications of tensor net- 
works for condensed matter problems. The relevant point to note is the two different categories 
of exact contractions vs. approximate solutions. The reduction we have just discussed is an 
exact one. That is, given ax = 2 tensor network discribing the unique solution to a 3SAT 
problem, there is a reduction problem to find the product state that delivers the same state. 
On the other hand, in condensed matter problems it is interesting to know whether a tensor 
network made of tensors of rank x can be faitfully approximated by a similar tensor network 
of lower rank x' ■ This is not an exact contraction and depends on the targeted accuracy 
of the tensor representation. In general the larger x is, the better the approximation be- 
comes. A quantitative analysis and characterization of scaling properties of a tensor network 
approximation as a funtion of x can be found in Ref. [241 125j . Further research devoted 
to determining the complexity of reducing the rank of a tensor network for special cases 
with a lattice geometry as the case of Ising of Potts models, are known to be NP-complete 
[261113113. 
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5 Conclusions 

Tensor networks provide a set of techniques to study by classical means interesting quantum 
states. These include universal resources for quantum computation, which can be expressed 
as tensor networks with x rank. Using a simple construction that encodes the solution of 
a 3SAT problem in a tensor network, \ip(Q,C)), we have shown that contracting arbitrary 
tensor networks, however loosely connected, may require exponential resources. From our 
construction, the number of solutions to an instance of the 3SAT problem can be directly 
read from the normalization (ip(Q , C)\ip(Q , C)) , hence the difficulty of this task is equivalent 
to a contraction of the network. The tensor network we have constructed is the ground state 
of a 3-local Hamiltonian. The hardness of finding its ground states implies the difficulty of 
cooling such a spin system. 

Let us emphasize that the difficulty of finding a solution to a 3SAT instance is not related 
to the construction of a tensor network for the state that encodes it. There is no need to resort 
to complicate minimization techniques. Rather, the difficulty is attached to the contraction of 
the tensor network. The problem of counting solutions and of finding explicitly one of them 
are solved with a polynomial number of contractions of the tensor network. It tourns out 
that a single contraction of the network counts the number of solution to the problem, thus 
solving a #P problem. The FNP problem of providing a single solution can be solved with a 
controlled sequence of contractions. Moreover, since the exact solution can be encoded in a 
X = 1 tensor network, finding a reduction in rank for the tensor network is also FNP. 

Since other hard problems reduce to the 3SAT problem addressed here, this result can 
be extended easily to an explicit construction of solutions to other hard problems like e.g. 
Exact Cover. Our tensor network encodes the relation between clauses in a set of tensors, 
resembling the message passing structure of other approaches used to solve huge instances 
of the 3SAT problem [11]. This provides an important link between techniques of both 
disciplines, suggesting the use of schemes developed in the study of 3SAT to improve the 
performance of contracting arbitrary tensor networks in Quantum Mechanics. 

Let us make a further link between the above result on computational complexity classes 
for the contraction of tensor networks and the possibility of constructing such tensor networks 
in the lab. In Ref. [28] a practical protocol was presented to create Matrix Product States in 
the laboratory. This idea has been further pursued in Ref. [29] in the context of PEPS. Our 
results show that the solution to 3SAT can be reduced to create the tensor network \ip(Q,C)) 
we have presented plus a set of n measurements. The solution to the problem would be read 
from the collapse of the state to one of the solutions which are superposed in \il>(Q, C)). As a 
consequence, it is expected that creating the quantum state representing the tensor network 
we have presented would require an exponentially long time. This shows that only tensor 
networks that capture some special symmetry in the geometrical connectivity of the problem 
might be amenable to experimental realization. 
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